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Many important phenomena occurring in nature can 
be described as growth processes. For such a reason 
these are attracting an increasing attention from theo- 
rists and experimentalists alike. It is well known [1-3] 
that simple models can lead to very rich dynamical be- 
haviours such as steady-state patterns, growing struc- 
tures, moving fronts, dynamical scaling, etc.. A great 
deal of study has been dedicated to the relaxation dy- 
namics of model A and B, in the Hohenberg and Halperin 
classification [1], where a single field evolves towards equi- 
librium with non conserved order parameter dynamics 
(NCOP) or conserved order parameter dynamics (COP), 
respectively. On the contrary, the Phase field model [4,5] 
which describes the coupling of a NCOP system with a 
diffusive field such as temperature or concentration seems 
not to have been completely explored, in spite of the fact 
that it displays a variety of interesting features. For ex- 
ample, one could mention the formation of solid-liquid 
fronts during the solidification process, a phenomenon 
widespread in nature and the possibility of observing a 
vast class of patterns according to the initial conditions 
and external fields. An advancing solid front releases 
heat which in turn diffuses and promotes the solidifica- 
tion of the undercooled melt. One also observes that the 
shape of the moving front is subject to an instability due 
to the coupling with the diffusive field whose result is 
the formation of tips and dendritic growth. Such a phe- 
nomenon is due to the presence of the diffusive heat-mode 
at the interface which conspires to destabilize the solid- 
melt boundary as predicted by Mullins and Sekerka [6]. 
In the present Letter we study theoretically the evolution 
of an ^-component version of the Phase field model and 
show that the coupling between the massless transverse 
modes and the diffusive field produce an instability which 
is the analogue of the Mullins-Sekerka instability. 



We shall consider a suitable generalization of the 
Hamiltonian of Refs. [4,7,8] to two coupled N compo- 
nent vector fields (j) = ((l)i(x,t), (I)ff(x,t)) and U = 
(Ui(x,t),...,UN(x,t)): 



H[cl>{x), U{x,t)] = J dx[^^{Vcl>y 



(1) 



where r and g are the standard quadratic and quartic 
couplings of the Ginzburg-Landau model and the last 
term represent a bilinear coupling between the field (j) 
and U . Stability imposes g > and w > 0. As far as the 
equilibrium properties are involved, one sees that the U 
field can be traced out, and the Hamiltonian (1) reduces 
to the familiar 0(N) spherical model [9]. The evolu- 
tion towards equilibrium is described by the two coupled 
Langevin equations: 



d(t)a(x,t) 

dt 



dUa{x,t) 

di 



6(t)a(x,t) 



Fr/V^ 



H[<j>,U] + ^„{x,t) (2) 



H[cj>,U]+Ux,t) (3) 



Equation (2) represents the evolution of a NCOP 4> cou- 
pled to a COP U , described by eq. (3). The noises 
appearing on the right hand sides of eqs. (2)-(3) have 
zero averages and are non conserved and conserved, re- 
spectively: 

{r,„{x,t)r,p{x\t)) = 2TjT^8„^p8{x - x')6(t-t') (4) 



{^aix, tMx',t)) = -2TfTu6a,pVH{x - x')6{t - t') 

(5) 

where Tj is the temperature of the final equilibrium state 
whereas Vjj and T^ are the kinetic coefficients. The 
noises ^„ and ry„ have vanishing mutual correlations. 

To make contact with physically motivated models we 
notice that in the scalar case N = 1 eqs. (2)-(3) reduce 
to the well known Phase field model, devised by Langer 
[2,5], to understand the solidification process starting 
from an undercooled melt. To this purpose one performs 
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the transformation U = u — [K^jw and identifies with 
the solid-liquid order parameter and u with the temper- 
ature field shift from the melting value and obtains the 
familiar order parameter and heat diffusion equations. 

In order to study the properties of the system it is 
sufficient to consider the three equal time real space 
connected correlation functions C^^(r,t) = [4>a{R + 
r,t)(^„{R,t)), C^u(r,t) = {<i>a(R + r,t)Ua(R,t)) and 
Cuu(f,t) = {Ua(R + r,t) Ua(R,t)) which are indepen- 
dent of the index a by symmetry. 

Within the N ^ co limit the dynamical equations for 
the correlations can be written in a closed form and the 
exact asymptotic behaviour explicitly extracted. In this 
limit their Fourier transforms, the structure functions, 
Cij(k,t), evolve according to the set of coupled equa- 
tions: 



i5CV^(M) 

2 dt 



dC^u(k,t) 
dt 



+ M^u{k,t)C^u{k,t) + T^Tf (6) 

= Mu4k,t)C^4k,t) 

+ (Muu(k,t) + M^^{k,t))C^u{k,t) 

+ M^uik,t)Cuuik,t) (7) 



l dCuu(k,t) 
2 dt 



Mu4k,t)C^u{k,t) 



+ Muu{k,t)Cuu{k,t) + TuTfk\ (8) 
The matrix elements M are given by 
M^^(k,t)= -T^k^ + r + gS), 

M^u(k,t)= -T^fi (9) 
Mu4,(k,t)= -Tufik\ Muu(k,t) = -Tuwk^ 

where S{t) = J^j^^^j^d'^k C^^{k ,t) contains a phe- 

nomenological momentum cutoff and d is the space di- 
mensionality. The model for r < fi'^/w, g > and 
c? > 3 displays a high temperature paramagnetic phase 
and a low temperature ordered phase, for Tf < Tg, the 
critical temperature of the model. As found in mod- 
els A and B, = (ji^/w - r)(d - 2)/(gA'^-^Kd) with 
1/A'd = 2''-'^TT''/^T(d/2) where T(x) is the usual Gamma 
function, and the structure factor reads 



hm C^^ik,t) = — ■ — —. 

t^oo '^'^ k^ + r + gb(oo) — /w 



(10) 



In other words at equilibrium the diffusive field produces 
a trivial renormalization of the bare parameter r. The 
quantity (r + gS(t) — fi'^/w) vanishes asymptotically and 
the small k structure function diverges. Dynamically, in- 
stead, the presence of U induces non trivial effects on 
the field (j) because it acts on a time scale longer than 



the noise field, characterized by a short correlation time. 
In the following we shall treat the two fields on equal 
footing. The line Tclji) separates a high temperature 
phase from a low temperature phase characterized by 
infinite correlation length and divergent fiuctuations at 
zero wavevector, due to the presence of massless Gold- 
stone modes. 

Since the behavior at Tf = is representative of the 
entire ordered phase at Tf < Tg we can neglect the noise 
term in the following analysis [10]. 

In the rest of the Letter we shall discuss the analyt- 
ical predictions on the dynamical behaviour. These are 
compared with the numerical solution of eqs. (6)-(8) for 
quenches at Tf = starting from uncorrelated random 
initial conditions, c? = 3, A = 1 and r < 0. 

In order to obtain a physical understanding of the re- 
laxation process we notice that the mode coupling term 
S(t) varies slowly with respect to the independent vari- 
able t. Thus a local linear analysis of the system of 
equations (2)-(3) for Tf = is sufficient to extract the 
three kinds of regimes which characterize the approach 
to equilibrium. The linearized system has, for each value 
of k, two linearly independent solutions, proportional 
to exp[w_|_(A;)t] and exp[w_(A;)t], associated with the two 
eigenvalues uj^(k) and ui-(k) of the matrix M. 

The larger eigenvalue uj^(k) controls the dynamics and 
a small k expansion is appropriate to extract the qualita- 
tive behaviour of the solutions. For temperatures below 
Tc and in the early regime S(t) is small. Thus, since 
(r -|- gS(t)) < one finds: 



j+(k) = T^\r + gS\- r^-r 



-gS] 



k^ + C4k^ 



(11) 



where C4 is a negative coefficient. One observes C^^(k,t) 
to display a peak around k = 0, growing in time with a 
power law ~ f^/"^, as in the pure NCOP, i.e. the longest 
wavelength fiuctuations grow faster. This is due to the 
fact that C^^(k,t) grows exponentially with a rate ap- 
proximately — 2r^(r -|- gS + k^). The coupling fi to the 
diffusive field has little effect in this regime and the sys- 
tem relaxes as if it had to reach an equilibrium char- 
acterized by r -|- gS = 0. However, this equilibrium is 
never achieved because of the presence of the field U . In 
fact, after a characteristic time Tf the value of r -|- gS(t) 
vanishes due to the build-up of correlations and the dis- 
persion relation becomes: 



ui+(k) = ^T^Tu/J-k 



TrjW 



k'^ + 0(k^). (12) 



Within this second regime the system evolves as if r = 
g = and we observe the growth of correlations with fi- 
nite wavevectors. This fact has remarkable repercussions 
on the structure factor, because we observe the appear- 
ance of a fastest growing mode with k = kf > 0. At 
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t = Tf , kf jumps discontinuously from to a finite value. 
The position km(t) of the peak of the structure factor 
C^^(k,t) then moves towards larger values of k. This 
regime describes formation of spatially periodic struc- 
tures. This mechanism reminds the Mullins-Sekerka in- 
stability where the solid-melt interface becomes unstable 
to finite wavelength fiuctuations and one observes pat- 
tern formation or dendrites. 

Since S(t) grows monotonically when r + gS(t) > the 
dispersion relation takes the form: 



r + gS 



c^k^ 



(13) 



with C4 < 0. Notice that the long wavelength modes 
become unstable because of the positive sign of the co- 
efficient of order k^ . In fig. 1 we display the dispersion 
relation uj^(k) for two typical times, t < tj and t > tj . 
Here there exist values of k for which LO+{k) is positive, 
and a small perturbation with wavevectors falling within 
this range will be amplified. From the quasilinear analy- 
sis one can extract the value at which the fastest growing 
mode appears for the first time. For small jj, this occurs 
at the wave vector: 



(14) 



From scaling we argue that the criterion for the insta- 
bility to occur is that the domain size L{t) ~ ^1/2 j^ygt 
exceed the characteristic length \j given by 2'irkJ^ , there- 
fore Tf ~ for small fi. In fig. 2 we display the numer- 
ical results for Tf vs. for small fi, the agreement is 
remarkable. Beyond the linear analysis one finds that the 
quartic coupling hinders the unbounded growth of fiuc- 
tuations. The non linear term, in fact, has a twofold role 
since initially it promotes the formation of the peak of 
the structure factor at finite wavevectors, but finally pre- 
vents it from growing indefinitely and eventually drives 
it towards zero wavevectors. 

At very late times, in the regime of the third kind, 
the peak finally tends towards k = and the width 
shrinks as the system tends to equilibrate. As t ^ co 
and r -|- gS — fi'^/w 0~, the coefficient of order k^ 
in eq. (13) vanishes and the system tends towards the 
COP static fixed point, displaying the typical behavior 
of conserved dynamics. Within this regime the domain 
size grows as L(t) ~ t^/^ and km{t) ~ {hi{t)/tYl^ as pre- 
dicted by the exact solution of the corresponding COP 
model without diffusion {jj, = Q). These predictions agree 
with the numerical solution. In fig. 3 the structure fac- 
tors at different times are shown. 

These findings are consistent with the existence of mul- 
tiscaling in the very late regime {t ^ Tf ), a particular 
form of scaling where the exponent is scale dependent. 
In particular the dynamical structure factor clearly dis- 
plays a multiscaling behavior as shown in C^^(k,t) ~ 



[i.(t)2A;„(t)2-<i]'/'(«^/«^™) [9] as shown in fig. 4. Such a 
multiscaling behavior follows from the competition of two 
marginally distinct lengths, namely the domain size L(t) 
and kj^. 

To conclude we have studied a model which provides 
an interesting new scenario. The observed instability at 
finite wavevectors is due to the coupling between low en- 
ergy Nambu-Goldstone modes and the diffusive modes. 
We believe this is a paradigm for the Mullins-Sekerka 
type of phenomena where the soft modes are represented 
by the capillary wave spectrum associated with the solid- 
melt interface and the diffusive mode is the heat trans- 
port. These two fields concur to destabilize the solid- 
melt boundary in analogy with our findings. On physi- 
cal grounds, one expects this kind of instability to occur 
during phase separation, because small droplets can dis- 
sipate heat more efficiently and reach rapidly thermal 
equilibrium due to the Gibbs-Thomson effect, whereas 
larger droplets try to dissipate energy faster by creat- 
ing bulges thus increasing the curvature. As the system 
cools down upon reaching equilibrium the typical length 
of the bulges, Xf, increases and diverges together with 
the average domain-size L(t). 

Our study of the present model is non perturbative 
and reveals an unexpectedly rich dynamical behavior: it 
displays scaling behavior in the early regime (t < Tf) and 
multiscaling [9] in the late regime (t ^ Tf). Moreover it 
displays a transient characterized by an instability which 
reminds in many respects the one occuring in solidifica- 
tion kinetics. It remains to be explored in greater detail 
the role of fiuctuations although the noise is thought to 
be irrelevant below the critical temperature of the model. 
We thank M. Conti and L. Peliti for illuminating discus- 
sions. This work was partially supported by grants from 
INFM and INFN. 
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FIG. 1. Dispersion relation ui^ (k) as a function of k for dif- 
ferent regimes: early regime (upper curve), and late regime 
(lower curve). 
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FIG. 2. Scaling of r/ as function of fi. The broken line 
has slope —2. The squares are obtained from the numerical 
solution with F,^ = 1, F^ = 5, r = —0.5, g = 1, w = 0.05 and 
d = 3. 
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FIG. 3. Structure factor C^^(k,t) for /t = 0.01 at succes- 
sive times: a) t = 500, b) t = 1500, c) t = 3000, d) t = 5000, 
e) t = 6000, f) t = 9000, g) t = 80000. The parameters are 
the same as fig. 2. 




FIG. 4. Shape function F{x) = k'^{t) C^^{xkm,t) as a 
function of x = k/km for different times for fi = 0.01. The 
non perfect data collapse signals multiscaling. In the inset we 
report the numerical value of the multiscaling exponent (p{x) 
defined in the text obtained from the best fit of C^^(xkm,i) 
as a function of L{tYkm{tY~ ■ The parameters are the same 
as fig. 2. 
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